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Abstract 

The propagation of electromagnetic plane waves with negative phase velocity (NPV) 
is considered in Schwarzschild-(anti-)de Sitter spacetime. It is demonstrated that 
NPV propagation occurs in Schwarzschild-de Sitter spacetime at lower values of the 
cosmological constant than is the case for de Sitter spacetime. Furthermore, we report 
that neither is NPV propagation observed in Schwarzschild-anti-de Sitter spacetime, 
nor is it possible outside the event horizon of a Schwarzschild blackhole. 


1 Introduction 

This communication concerns the propagation of electromagnetic planewaves with negative 
phase velocity (NPV) in curved spacetime. The phase velocity is called negative if the time- 
averaged Poynting vector and the wavevector are oppositely aligned [1]. The phenomenon 
of negative refraction — which has lately generated considerable excitement in the elec¬ 
tromagnetics and materials research communities [2] — follows as a direct consequence of 
NPV propagation. The prospect of technological applications, such as lenses with extremely 
low distortion, has prompted intense efforts by experimentalists and theoreticians directed 
towards the development of NPV-supporting artificial metamaterials [3]. 

The importance of NPV propagation in astrophysical scenarios has recently been em¬ 
phasized. We have previously shown that vacuum, in association with certain spacetime 
metrics, can support NPV propagation [4, 5]. For example, the metric of de Sitter spacetime 
admits NPV propagation, whereas the anti-de Sitter metric does not [6]. Thus, an exper¬ 
imental means of distinguishing between anti-de Sitter and de Sitter spacetimes is offered 
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via NPV propagation. Furthermore, regions supporting NPV propagation are to be found 
within the ergosphere of a rotating black hole [7]. In this context, comparisons between NPV 
propagation and superradiant scattering are noteworthy. While both phenomenons involve 
negative energy densities, there are significant differences in terms of directional properties 
and frequency bounds [8]. 

In the present study we consider planewave propagation in Schwarzschild-(anti-)de Sitter 
spacetime. Our analysis is based on a formulation originally proposed by Tamm [9], wherein 
the formal analogy between electromagnetic propagation in gravitationally affected vacuum 
and propagation in a (fictitious) instantaneously responding medium is exploited [10] . 


2 Schwarzschild— (anti—)de Sitter spacetime 

The static Schwarzscliild-(anti -)de Sitter metric, with the signature (+, —, —, —), is conven¬ 
tionally expressed in spherical coordinates as [11, 12, 13, 14, 15] 

ds 2 = (1 — f)dt 2 — Y~~y d 1 " 2 ~ r2 (d@ 2 + sin 2 $ dqf 2 ) , (1) 

wherein / = (2 GMr^ 1 + Ar 2 /3)/c 2 , M is the mass of the black hole, c is the speed of light 
in vacuum in the absence of a gravitational held, G is the gravitational constant, and A 
is the cosmological constant. If the cosmological constant A is positive, the spacetime is 
called Schwarzschild-de Sitter spacetime, whereas Schwarzschild-anti-de Sitter spacetime is 
characterized by A < 0 [16]. Let us note that (1) reduces to 

(a) the Schwarzschild metric when A = 0, and 

(b) the cosmological metric (de Sitter or anti-de Sitter, according to whether A is positive 
or negative, respectively) when M = 0. 

In terms of Cartesian coordinates x = rsin#cos0, y = r sin#sin0, and z = r cos#, the 
metric (1) is represented by g a g as 4 


/ 1-/ 


[ 9a/3 ] ~ 


V 


o 


-l - 


o 

fx 2 

r 2 (1 - /) 
fxy 


0 

fxy 


r (1 -/) 


r (1 -/) 


fxz 


r (1 -/) 


fv 2 

r 2 (1 - /) 
fyz 


o 

.fxz 


r (1 - /) 


fyz 


r (1 - /) 


r (1 - /) 


fz 2 


\ 


_ r 2 (1 - /) / 

4 Roman indexes take the values 1, 2 and 3; while Greek indexes take the values 0, 1, 2, and 3. 
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and its inverse g af3 as 
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Following common practice [10, 17, 18, 19] , the electromagnetic response of vacuum in 
curved spacetime may be described by the constitutive relations of an equivalent, instanta¬ 
neously responding, medium as per 


D = £o2'E 1 

B = g7l-H j’ 


wherein SI units are implemented. Here, e 0 = 8.854 x 10 12 F m 1 , /i 0 = 4 -k x 10 12 H m 
and 7 is the 3x3 dyadic equivalent of the metric [ 7 ^] with components 


0 ^ab 



( 5 ) 


3 Piecewise uniform approximation 


A global description of Schwarzschild-(anti-)de Sitter spacetime is provided by the consti¬ 
tutive relations (4). Let us partition the global spacetime into adjoining neighbourhoods. 
At an arbitrary location (x,y,z), we consider the neighbourhood 77 which is taken to be 
sufficiently small that the nonuniform metric 7 a b may be reasonably approximated by the 
uniform metric 7 ^ [7]. By stitching together solutions from adjoining neighbourhoods, we 
formulate the global solution. This piecewise uniform approximation technique is commonly 
employed in solving differential equations with nonhomogeneous coefficients [20]. Thus, we 
have the uniform 3x3 dyadic representation 


7 = [7 ab 
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at 1Z , with 


(7) 
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and r 2 = a ; 2 + y 2 + 5 2 . We note that det 7 = ^1 — / 

4 Plane waves in 

We seek planewave solutions 

E = Re {E _ 0 exp [i(k’r — tot )}} , H_ = Re { H_ 0 exp [i (k'r — tot )]} , ( 8 ) 

to the source-free Maxwell curl postulates 

r\ r\ 

X7xE + ^-B = 0, VxU-^ = 0 (9) 

— dt~ ~ — dt~ ~ y J 

in 1Z. Here, the wavevector is denoted by k and the position vector within the neighbour¬ 
hood containing (x, y, z ) is represented by r; to and t are the angular frequency and time, 
respectively. Let us emphasize here that r and t are independent of TZ. 

The amplitudes E 0 and H _ 0 are complex-valued with i = I. Upon combining ( 8 ) and 
(9), we find after some algebraic manipulation that 

W • E_ 0 = 0 , (10) 

where 

W = ^kgdet 7 — k • 7 • kj 1+ kk • 7 , (11) 

and k 0 = to^/e 0 y 0 . Thus, the dispersion relation det [ W ] = 0 emerges, which can be recast 
as 

/c 2 det 7 (kl det 7 — k • 7 • k j = 0 . ( 12 ) 

Clearly, the wavevectors must satisfy the condition 

k'~f'k = kl det 7 , (13) 

if 7 is nonsingular. 

Let us consider now the eigensolutions of (10). In view of (13), we have 

kk • 7 • E_ 0 = 0 ; (14) 

hence, it follows that E_ 0 is orthogonal to k • 7 . The spacetime metric described by (1) is 
spherically symmetric. Therefore, there is no loss of generality in choosing the wavevector 

k = ku z , (15) 
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with u z being the unit vector lying along the z Cartesian axis. Thereby, 

k*Z = + 72 u y + 73 u z ), (16) 

where 

~ fxz _ /P _ /* 

7l - —7-7\—r > 72 - —7-Z\—r , 73 - -— , (17) 
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and and u y being unit vectors lying along the x and y Cartesian axes, respectively. 
The two linearly independent eigenvectors 

£i = l2lL x - 7i u y , £2 = 7i73^x + 727s u y - (7 2 + u z 

satisfy (14); hence, we have the general solution 

E 0 = C'ie 1 + C2e 2 , 


(19) 


( 20 ) 


with Ci and C' 2 being arbitrary complex-valued constants. The corresponding expression 
for H_ 0 follows from the Maxwell postulates as 
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To calculate the wavenumbers we turn to the dispersion equation (13). For k aligned 
with u z we obtain the /e-quadratic expression 
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from which the wavenumbers 


k = ±A)o 
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straightforwardly emerge. 

The requirement that fcel imposes the condition 


1-./ ) fz > o; 


(24) 


in other words, both (1 — /) and f z must have the same signs for propagating planewave 
solutions. 
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5 NPV Condition 

The propagation of planewaves with NPV is signalled by the inequality [5, 7] 


k'{P)t< 0, (25) 

where (P_)t — (1/2) Re {E_ 0 x H_* 0 } is the the time-averaged Poynting vector. The general 
solution (20) and (21) delivers 


(£>< = 7^2— I ic-ii 2 (i -/) + |C 2 | 2 /: 


l 1 *' V V ■ "J - 1 X - 2 ' 

By virtue of (19), the orientation of ( P_) t is provided by the vector 

Li x e 2 = ( 7 ? + 7i) (fii* + 72« y + 7s.«*) • 

Thus, 
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Hence, by exploiting (24), we see that NPV arises as a consequence of f z < 0 ; i.e., NPV 
propagation occurs provided that 


6 GM 3c 2 
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(30) 


Let us emphasize that while (30) has been derived for the neighbourhood 77, the location of 77 
is arbitrary within Schwarzschild-(anti-)de Sitter spacetime. Therefore, the NPV inequality 
(30) applies generally. 

In order to contextualize the NPV condition (30), it is helpful to introduce the event 
horizon for a Schwarzschild black hole which lies at r = r Sch = 2 GM/c 2 and the event 
horizon for de Sitter spacetime which lies at r = r deS = c^/zJX [21]. Thereby, we see that 
(30) yields the following sufficient conditions for NPV propagation: 


' > r des — r sch) for A > 0 
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6 Concluding remarks 


The general condition (30) has been derived for NPV propagation in Schwarzschild-(anti-)de 
Sitter spacetime. Analysis of this condition allows us to make the following conclusions: 


(a) If A = 0 (i.e., the spacetime is described by the Schwarzschild metric), then NPV 
propagation is indicated by the inequality 


f sa 

r > z 2 ’ 


(32) 


Therefore, we deduce that NPV propagation is not possible outside the event horizon 
of a Schwarzschild black hole (i.e., for r > r Sch ). 


(b) If M — 0 (i.e., the spacetime is (anti-)de Sitter), then NPV is signalled by the inequality 

3r 2 

A>—, (33) 

z z 

in agreement with earlier results [6]; i.e., anti-de Sitter spacetime does not support 
NPV propagation. Furthermore, we observe that it is the region of spacetime which 
lies outside the de Sitter event horizon which supports NPV propagation. 

Let us analyze the NPV inequality (33) in view of the uniform approximation imple¬ 
mented in considering the neighbourhood 1Z. The linear dimensions 5 of 7 Z are taken to 
be small relative to the global spacetime curvature and large compared to the electro¬ 
magnetic wavelengths as given by 2ir/k. Since the Ricci scalar R — which provides a 
measure of the inverse radius of spacetime curvature squared — is given by R = 4A/c 2 
for de Sitter spacetime [22], we have 
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with p being a proportionality constant, 
requires 

|A| « 


Hence, the partition of global spacetime 
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The inequalities (33) and (34) are mutually compatible as the linear dimensions 5 of the 
neighbourhood 1Z are chosen independently of the 5 coordinate specifying the location 
of n. 


(c) By comparing (30) and (33), it is clear that in Schwarzschild-de Sitter spacetime NPV 
propagation occurs for smaller (positive) values of A than is the case for de Sitter 
spacetime. 

(d) Our numerical investigations have shown that the NPV inequality (30) is not satisfied 
for A < 0 (i.e., Schwarzschild-anti-de Sitter spacetime) in physically probe-able regions 
of spacetime. 
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